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We investigate numerically the resistive transition in the two-dimensional XY gauge glass model. 
The resistively-shunted junction dynamics subject to the fluctuating twist boundary condition is used 
and the linear resistances in the absence of an external current at various system sizes are computed. 
Through the use of the standard finite-size scaling method, the finite temperature resistive transition 
is found at fcsTc — 0.22(2) (in units of the Josephson coupling strength) with dynamic critical 
exponent z — 2.0(1) and the static exponent = 1.2(2), in contrast to widely believed expectation 
of the zero-temperature transition. Comparisons with existing experiments and simulations are also 
made. 

PACS numbers: 75.10.Nr, 64.70.Pf, 74.50. -hr, 74.25.Fy 



Since the prediction Qfr-the vortex glass (VG) phase 
in the high-Tg materialspcl the properties and the exis- 
tence of. iiCI phase have drawn intensive interest for a 
decade.B"t3 In experimentaJLl the VG phase has a sig- 
nificant practical importance since the high- Tc materials 
can be truly superconducting in this phase where vor- 
tices are completely frozen at random positions, while 
the Abrikosov vortex lattice in pure type-II superconduc- 
tors dissipate energy at any amount of external currents. 
To study VG from a theoretical point of view simplified 
discretized models such as the XY gauge glass and the 
XY spin glass have been widely used. For the gauge 
glass model, there is a growing consensus that the low 
critical dimension is less than 3, which means that there 
exists VG phase at finite temperatures in three dimen- 
sions (3D). This has been confirmed by ihe. defect wall 
energy calculations at zero tempeSjaitareEl^El as well as 
by finite-temperature simulations. t3"llil Although the is- 
sue about the validity of using the gauge glass model to 
describe real bulk high-Tj, materials is not completely set- 
tled (e.g., the former is isotropic with a vanishing net ex- 
ternal magnetic field, while the latter is anisotropic with 
nonzero net magnetic field), experiments on high- Tc ma- 
terials like YBColrla and BiSCCOffl also have yielded 
results in accordance with the existence of a VG phase 
at finite temperatures. 

In 2D, all existing defect wall energy calculations have 
unanimously rayealed that the stiffness exponent has a 
negative value ,BEilll3 which has been confirmed in Ref. |l^ 
through a simple argument but has been interpreted as 
an artifact of ubiquitous spin wave fluctuations. Finite- 
temperature simulations have so far yielded contradict- 
ing results: On the one hand, zero-temperature transi- 
tion has been concluded imrn. Monte-Carlo (MG) sim- 
ulations in earlier studiesllj'E£lEZl and from the current- 
voltage {IV) characteristics, in resistively-shunted junc- 
tion (RSJ) simulations .ESEj On the other hand, there 
have been evidences of thCpSjiite-temperature transition 
from the IV characteristic's and from the divergence of 
the relaxation time scaleEil in RSJ dynamic simulations, 
and from the very recent extensive finite-temperature 



MG simulations. Also, Langevin dynamic simulation 
of a vortex model with quenched impurities, which is 
closely related with the gauge glass model— also has 
found the finite-temperature transition in 2D.E3 To make 
things more complicated, absence of finite-temperature 
VG phase has been reported in early experiments on 
very thin YBCO filmsp whereas a recent study on thin 
YBCO films, where the coiselation length in c-direction 
exceeds the film thickness J3 as well as experiments on 
highly isotropic BiSGGO [Bi(2:2:2:3)]l!l has obtained re- 
sults which can be interpreted as indications of finite-T 
VG phase in 2D. Gonsequently, we believe that the ques- 
tion about the existence of the finite-T VG transition is 
not resolved yet. To our knowledge, a detailed numerical 
study of transport properties of 2D gauge glass, pipdel 
has not been performed: Existing RSJ studiesEa'Ej'LiJ did 
not consider the-Jmitc-size scaling in a proper way and 
in some of therr£3't2l the temperature range used did not 
cover the expected transition temperature (Tc « 0.22 in 
Refs. |l6| and |l[ and T^ « 0.15 in Ref. |o|). 

In this paper, we use RSJ dynamic simulation to inves- 
tigate the resistive transition in the 2D XY gauge glass 
model in the absence of external currents. Since the volt- 
age is always nonzero in the presence of a finite current 
because of the nucleation process, vanishing of the linear 
resistance, which is defined in the limit of zero external 
current, is the appropriate definition of a superconduc- 
tor. The fluctuating twist boundary [eondition (FTBG) 
for Langevin-type dynamic simulationcJ made it possible 
to calculate the linear resistance with preserved period- 
icity of the phase variables (similar methods have also 
been used in Ref. |2|). We here use FTBG to calculate 
the linear resistances, and then, through the standard 
finite-size scaling analysis (see, e.g., Refs. || and ^), a 
resistive transition at a finite temperature is concluded. 

The Hamiltonian of the 2D L x L^JCY. gauge glass 
model subject to FTBG is written asoEZrEa 

H^-jJ2 cos{9, - 6, - r,y • A - A„ ), (1) 
where J is the Josephson coupling strength, the summa- 
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tion is over all nearest neighboring bonds, and 9i (with 
periodicity 9i = Oi+Lit = ^i+Ly) is the phase of the su- 
perconducting order parameter at site i. The twist vari- 
able A — (A^;, Ay) measures the global twist of phase 
variables in each direction, and {— x, y) is the unit 
vector from site i to the nearest-neighbor site j (the lat- 
tice spacing is taken to be unity). The XY gauge glass 
model is characterized by the quenched random variable 
Aij uniformly distributed in [— 7r,7r), whereas in the XY 
spin glass model Aij can have binary values and tt with 
equal probability. Recently XY spin glass model draws 
attention in relation to the 7r-junctions due to the d-wave 
symmetry of order parameters. 

The equations of motion for RSJ dynamics subject to 
FTBC are determined from the local and the global cur- 
rent conservations and are written in dimcnsionless forms 
(see Refs. p3 and |23 for details): For phase variables 
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where is the 2D square lattice Green's function and 
the primed summation is over four nearest neighbors k 
of site j, while for the twist variables 
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The thermal noise terms satisfy (rjij) = (^^a) 
iVijil'^) = (f?A,?/A„) = and 
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where (■ • ■) is the thermal average, and the time t and 
the temperature T have been normalized in units of 
fi^ / Ae^R^J with the shunt resistance i?o, and J/kB, re- 
spectively. The equations of motion were integrated 
numerically by using the efficient second-order Runge- 
Kutta-Helfand-Greenside algorithrrO with the discrete 
time step At = 0.05. For each disorder realization we 
neglected an initial stage of the time evolution up to tg, 
and the measurements were made between and io + ^m- 
For L = 4, 6, 8, and 10, disorder averages were performed 
over Ns = 400, 200, 100, 50 samples, respectively. We 
verified that to = 10^ and = 5 x 10^ (corresponding to 
10^ time steps) were sufficiently large; it is more impor- 
tant to increase the number of sample averages to have 
more reliable results than increasing tm further. The 
highest temperature T = 0.30 in the present work cor- 
responds to the lowest temperature in Ref. |l8|, and this 
makes it necessary to use much larger value of i,„ than 
in Ref. ^ [by factor 0(10^)]. 

The linear resistance in x direction is given by the 
Nyquist formula:E3 
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with the disorder average [■ • ■] and the voltage across the 
whale sample Vx = —Li^^, which is then approximated 
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[([A,(e)-A,(0)]^ 
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for sufficiently large 0. For convenience, we define 

_ [([A.(e) - A.(0)]2)] + [([A,(e) - A,(0)]2) 
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and calculate the linear resistance R from the least- 
square fit to the form f{Q) — i?8. 

We plot in Fig. for L = 6 at T = 0.22 (a) the time 
evolution of iAy for a given disorder realization, and (b) 
/(O) obtained from the average over 200 samples. The 
full line in Fig. |l|(b) obtained from the fit determines the 
value of R. Figure ^ shows R obtained in this way as 
a function of L at T = 0.30, 0.25, 0.23, 0.22, 0.21, 0.20, 
and 0.15. For each L and T, we divided the total number 
of samples, Ns^ into five groups, each of which contains 
A^s/5 samples, and R was calculated for each group, lead- 
ing to the estimation of error bars (standard deviations) 
in Fig. ^ The linear resistance shows clear change of cur- 
vature around T « 0.22: In the high-temperature regime 
R appears to saturate as L is increased, while in the low- 
temperature regime R drops down more rapidly as L is 
increased. To see this behavior more clearly, we also dis- 
play in the inset of Fig. | R{T = 0. 15) /i?(T = 0.30) as 
a function of L in log scales, which does not exhibit any 
sign of saturation in terms of L: On the contrary, the 
downward curvature implies that in the thermodynamic 
limit i?(T = 0.15) vanishes. These observations lead to 
the conclusion that the 2D XY gauge glass model has 
the finite-temperature resistive transition. 

To obtain more precise value of Tc, we proceed to the 
finite-size scaling analysis of the linear resistance. In the 
dynamic scaling theory for the usual second-order phase 
transition in 2D,cl the linear resistance has the standard 
finite-size scaling form: 



R{L, T) = L-^p f (T - Tc)L^''' 
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where p{x) is the scaling function with the scaling vari- 
able X, z is the dynamic critical exponent, and the static 
critical exponent v is defined by the divergence of the cor- 
relation length, i.e., $ ~ \T -T^^ . MT = Tc, R^ L"^ 
[see Eq. (|^)] and from Fig. ^we obtain z « 2.0. Although 
this value is much smaller than the values usually mea- 
sured in experiments (z = 4 ~ 9), it should be noted 
that the similar values have also been observed in vari- 
ous numerical simulations with RSJ dynamics and MC 
dynamics of 2D gauge glass model {z « 2.2 in Ref. pO|and 
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z w 2.4 in Ref. |6[ respectively) as well as the Langevin- 
type relaxational dynamics of vortices in the presence of 
quenched impurities {z « 2.1, Ref. 23). Furthermore, a 
recent analytic calculation based on the dynamic renor- 
malization group method also has found z — 2 for 2D 
gauge glass model vifith a purely relaxational equations of 
motionE3 In contrast, widely believed zero-temperature 
transition in 2D gauge glass model implies that the linear 
resistance should have Arrhenius form and thus corre- 
sponds to z = oo. We aisa,tried to fit our data in Fig. ^ 
to the Arrhenius formllallj (see Fig. H) and found that 
the thermal activation barrier, which is proportional to 
the slopes in Fig. ||, strongly depends on the system size. 
Furthermore, R is found to deviate from the Arrhenius 
form in a systematic way at the lowest temperature at 
all system sizes. 

Figure § shows the scaling plot Eq. (0) with z — 2.0, 
Tc = 0.22, and z/ = 1.2. We tried to vary parameter 
values used in the scaling plot and concluded 



T, = 0.22(2), 
z = 2.0(1), 
= 1.2(2), 
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where numbers in parenthesis denote errors in the last 
digits, and Tc ~ 0.22 is in a good agreement with Refs. |l^ 
and p^. It is worth mentioning that the same method 
(calculation of the linear resistance with FTBC accom- 
panied by the finite-size scaling analysis) has leadj-tp a 
very precise determination of Tc in 3D XY model.Ea In 
recent studies of 2D iLY_pauge glass model based on the 
same RSJ dynamics Jla^Ej the finite-size scaling analysis 
has not been used and the temperature range did not 
cover the critical temperature (Tc « 0.22), leading to the 
different conclusion of the zero-temperature transition. 

Recently, the possibility of the quasi-long-rauae glass 
order with the vanishing glass order parameterE3 in 2D 
gauge glass model has been suggested in Ref. where it 
has been argued that the correlation length ^ diverges in 
the whole low-temperature phase. Although our scaling 
plot in Fig. ^was obtained on the assumption of the usual 
second order transition, where ^ is finite both below and 
above Tc, it is still plausible to have the quasi-long-range 
glass order: In this case, the quasi-criticality should be 
reflected in the temperature-dependent dynamic critical 
exponent below Tc with R{L,T) ~ L'^C^', while the 
high-temperature phase still obeys the scaling form in 
Eq. (|^). Our current results cannot rule out this possi- 
bility, and if we neglect the smallest size i = 4 in Fig. ||, 
R in the low-temperature phase indeed appears to ex- 
hibit the simple algebraic form R with increasing 
z as T is decreased, in accord with the idea of the quasi- 
long-range glass order. 

In conclusion, we studied numerically the resistive 
transition in the 2D XY gauge glass model by using 
the RSJ dynamic equations subject to the fluctuating 
twist boundary condition. The standard finite-size scal- 
ing analysis applied to the linear resistances lead to the 



strong evidence of the finite-temperature resistive transi- 
tion at Tc = 0.22(2) with the dynamic critical exponent 
z = 2.0(1) and the static exponent v = 1.2(2). However, 
the nature of this finite-temperature resistive transition 
needs to be investigated more in detail. 
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FIG. 1. (a) Time evolution of LAy{t) for a given disor- 
der realization for L = 6 at T = 0.22. (b) /(0) versus O 
[see Eq. 1^] obtained from 200 sample averages for L — 6 
a,t T — 0.22 (open circles). Least-square fit to the form 
/(6) = RO (full line) determines the linear resistance. 




FIG. 3. Arrhenius plot of the linear resistance: Ini? versus 
1/T. R shows a systematic deviation at the lowest tempera- 
ture for all system sizes used. 
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FIG. 4. Finite-size scaling plot of the linear resistance: 
RL'' versus the scaling variable (T - Tc)L^''" with Tc = 0.22, 
FIG. 2. Linear resistance R versus system size L at tern- = 1.2, and z = 2.0. (All data points in Fig. || have been 

peratures T = 0.30, 0.25, 0.23, 0.22, 0.21, 0.20, and 0.15 used.) 
(from top to bottom). It is shown that the behavior in 
the high-temperature regime {T > 0.25) is different from 
that in the low-temperature regime (T < 0.15), imply- 
ing the existence of the resistive transition near T — 0.22, 
where the power-law behavior R ~ also determines 

the dynamic critical exponent z ~ 2. Inset: The ratio 
R(T = 0.15)/i?(r = 0.30) versus L. It is exhibited that the 
ratio shows downward curvature as L is increased, which im- 
plies that R{T = 0.15) vanishes in the thermodynamic limit 
and that the resistive transition occurs at Tc > 0.15. 
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